Vortex dynamics in superconducting channels with periodic constrictions 
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Vortices confined to superconducting easy flow channels with periodic constrictions exhibit re- 
versible oscillations in the critical current at which vortices begin moving as the external magnetic 
field is varied. This commensurability scales with the channel shape and arrangement, although 
screening effects play an important role. For large magnetic fields, some of the vortices become 
pinned outside of the channels, leading to magnetic hysteresis in the critical current. Some channel 
configurations also exhibit a dynamical hysteresis in the flux-flow regime near the matching fields. 

PACS numbers: 74.25.Wx, 74.78. Na, 74.25. Sv 



Vortices flowing through nanofabricated easy flow 
channels in superconducting films provide a useful system 
for studying the dynamics of interacting particles mov- 
ing in tailored confining potentials. The general problem 
of interacting particles in confined geometries is impor- 
tant in a variety of physical systems, including colloids 
flowing through microchannels [1] and Wigner crystals 
in the presence of constrictions 0). With appropri- 
ate asymmetries, such tailored potentials can also form 
model systems for studying ratchet dynamics, with ap- 
plications ranging from superconducting devices to inves- 
tigations of biomolecular motors Q ■ The fabrication of 
weak-pinning channels for guiding vortices through su- 
erconducting films at the nanoscale is well established 
Such channels have been employed in a variety of 
investigations of vortex dynamics at relatively large mag- 
netic fields, typically greater than 10 3 Oe, including ex- 
periments on mode locking Q and melting in confined 
geometries [7|. 

Recent advances in nanofabrication have enabled im- 
plementations of artificial periodic vortex pinning lat- 
tices in superconducting films. These are typically pro- 
duced with arrays of either nanoscale holes through the 
film d, Q or magnetic dots underneath the film fioj j . 
Such structures result in a substantial magnetic field- 
dependence to the critical current, which is related to 
the threshold force required to cause vortex motion. The 
critical current typically exhibits commensurate behav- 
ior with maxima when the magnetic field corresponds to 
an integer number of vortices per pinning site. For fields 
away from these matching points, the dynamics of inter- 
stitial vortices, which are not located on the strong pin- 
ning sites but rather are more weakly confined through 
interactions with the strongly pinned vortices, lead to 
lower critical currents. A variety of experiments have 
been performed on such pinning arrays in recent years, in- 
cluding studies of the pin ning -strength dependence [111 ], 
quasiperiodic lattices [FJ, LL3J , and structures with ran- 
dom dilutions of pinning sites There have been 
many simulations of vortex dynamics in these periodic 
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FIG. 1: (Color online) (a) Scanning electron micrograph of 
three channels with periodic diamond constrictions, (b) Strip 
layout, along with channel and magnetic field orientation. 



pinning systems as well [130]. 

In this article, we describe measurements of vortex 
dynamics in weak-pinning channels that contain peri- 
odic constrictions at small magnetic fields, generally less 
than 10 Oe. Thus, this involves considerably smaller 
fields than much of the previous work on vortex matter 
in unstructured weak-pinning channels. The nature of 
the channels provides pathways for the easy flow of vor- 
tices, while the lattice of periodic constrictions results in 
strong matching effects with substantial enhancement of 
the critical current I c at certain values of the external 
magnetic field H a . Although we do not image our vortex 
distributions directly, we can determine that over much 
of the field-range of our measurements, all of the vortices 
are confined to the channels, with the dynamics deter- 
mined solely by the channel geometry, screening currents 
in the film, and interactions between vortices. Thus, in 
this field regime there is no distinction between pinned 
and interstitial vortices. At larger H a vortices can enter 
the regions outside of the channels where they become 
strongly pinned and do not participate in the flux-flow. 
Instead, these pinned vortices alter the potential for the 
vortices that are confined to the channels and lead to an 
irreversibility of I c (H a ). Besides the magnetic-field hys- 
teresis in I c for large H a , we often observe a dynamical 
hysteresis in the vicinity of the matching fields in the 
current- voltage characteristics (IVCs) themselves. 

Following the scheme in Refs. [5H7(, we fabricate our 
channels from bilayer films of a 200 nm-thick layer of 
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amorphous-NbGe, an extremely weak-pinning supercon- 
ductor (T c NbGe = 2.93 K), and a 50 nm-thick NbN layer, 
with relatively strong pinning (T c NbN w 10 K), on a Si 
substrate. The channels are defined with electron-beam 
lithography, followed by a reactive ion etching process 
to remove the NbN, resulting in weak-pinning channels 
for vortices to move through easily. The channels are ar- 
ranged across a 50 /im-wide strip, with H a oriented along 
the thin axis of the strip (Fig. [1]). The strip pattern con- 
tain pairs of probes for coupling to a room-temperature 
low-noise amplifier for sensing the voltage drop V along 
the strip due to vortex motion through the channels. A 
transport current driven through the strip with an ex- 
ternal supply generates a transverse Lorentz force on the 
vortices. Between each pair of voltage probes is an ar- 
ray of 15 identical channels with inter-channel spacing 
s. Each channel contains a periodic chain of cells de- 
fined by diamond-shaped constrictions, all of which are 
3.2 /im across at the widest point and 700 nm wide at 
the constriction, with a period along the channel p. We 
have measured sets of such channels with five different 
combinations of (s,p). 

We perform our measurements with the strip immersed 
in a pumped helium bath. Our results presented here 
were obtained at temperature T between 2.61 K and 
2.90 K (89% - 99% of T c NbGo ). We can apply the stan- 
dard dirty-limit expressions to estimate the relevant su- 
perconducting parameters of the a-NbGe and NbN films. 
For the a-NbGe, the coherence length £ varies between 
20 — 80 nm over the range of T, thus, the vortex core size 
is always much less than the smallest dimension of the 
channels and the vortex cores are essentially point-like. 
On the other hand, the penetration depth is quite large, 
and the thin-film screening length, Aj^ = 2A 2 /d, where d 
is the film thickness, ranges between 40 — 370 /im for the 
a-NbGe. In the NbN that forms the banks between the 
channels, w 8 /im with little temperature variation 

since T/T c NbN < 1 [18]. Thus, the circulating currents 
for a vortex in a NbGe channel extend along many, if not 
all, of the diamond cells in that particular channel and 
penetrate roughly 8 /jm into the NbN banks on either 
side of the channel. Because A^ b is much greater than 
both A^ bN and the width of the channels, vortices will 
be confined to the channels and the shape of the channel 
walls will play an important role in distorting the cir- 
culating currents around each vortex. As in our earlier 
measurements of ratchet dynamics with asymmetrically 
distorted weak-pinning channels 18], by controlling the 
channel wall shape, it is possible to tailor the confining 
potential for a vortex in the channel. 

We characterize the transition from the static state to 
a dynamical flux-flow regime by measuring the critical 
current I c in the conventional way, that is, by monitor- 
ing the current-voltage characteristic, then applying a 
1 /iV criterion. We drive the vortices with 200 cycles of 
a bias current sinusoid at 210 Hz, then average the re- 
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FIG. 2: (Color online) (a) Measurement of I c (H a ) for p — 
2 /Ltm, s = 20 /nm, T — 2.78 K for a complete field cycle as de- 
scribed in the text with arrows and legend indicating different 
portions of magnetic field sweep, (inset) I c {H a ) for 0.5 fim- 
wide uniform channels (s = 20 /im, T = 2.78 K). (b) I c (H a ) 
for same channel parameters as the main figure for different T 
as indicated, scaled by corresponding i c (0). (c) Correspond- 
ing Fourier transform magnitudes with vertical dotted line 
indicating location of uijj for this particular configuration of 
channels. 



suiting voltage response. We generate H a with a super- 
conducting Helmholtz coil and a /i-metal shield reduces 
the background magnetic field below 13 mG. For each 
measurement sequence, the strip was heated to ~ 17 K, 
well above T c of both the NbGe and NbN films, and was 
then cooled in H a = 0, while we subsequently increased 
H a at the measurement temperature. 

Measurements of the field dependence I c (H a ) yield in- 
formation about the vortex dynamics in the channels. 
For comparison, we fabricated a set of 0.5 /xm-wide uni- 
form channels, thus, with no constrictions, and measured 
I c {H a ) [Fig. dja) (inset)]. The response is similar to that 
characteristic of an edge barrier for a thin, weak-pinning 
superconducting strip in a perpendicular magnetic field, 
where the entry of vortices at the strip edge is deter- 
mined by the distortion of the current density across the 
width of the strip lj| 20j]. For a standard edge barrier, 
I c {H a ) follows two different regimes: for H a near zero, I c 



decreases linearly with H a , when vortices enter the strip 
at one edge and are immediately swept across the entire 
strip width; for larger H a , I c oc H^ 1 , where the external 
magnetic field is large enough to push vortices into the 
strip, even for transport currents less than I c . 
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FIG. 3: (Color online) (a) I c {H a ) curves for strips with dif- 
ferent p for s = 20 fim, T = 2.78 K. (b) p-dependence of w% ; 
dashed line is a guide to the eye (slope = 0.46 Oe _1 /im _1 ); 
dotted line has slope 2A 1 1 bN /<I>o (= 0.77 Oe - 1 /inT 1 ). (c) 
Ic{H a ) curves for different s, with p = 2 /im, T = 2.78 K. (d) 
s-dependence of uj° h \ horizontal dashed line at 2pA^ bN /<I>o, 
dotted line has slope p/&o- 

The presence of diamond-shaped constrictions in the 
channels results in pronounced oscillations in I c (H a ) on 
top of the edge barrier response [Fig.^a)]. For this mea- 
surement, H a was increased from to 6.2 Oe, then re- 
duced through to —6.2 Oe, and finally returned to 0. 
The complete reversibility of I c (H a ) for this field-cycle 
indicates that all of the vortices are confined to the chan- 
nels, as one would expect a reversible I c (H a ) for a pure 
edge barrier. In contrast, if vortices had entered the 
strong-pinning NbN, one would expect to observe hys- 
teresis in I c (H a ). The oscillations in I c (H a ) can be ob- 
served over a wide range of T [Fig.^b)], with the relative 
height of the peaks increasing as 

jy T NbGe approaches 
1. A Fourier transform of the I c (H a ) data [Fig. dfc)] 
shows that the characteristic frequency of these oscilla- 
tions, uP H = l/AH a , as identified by the vertical dotted 
line in the figure, is independent of T in this range, indi- 
cating that the commensurability is determined primarily 
by the channel geometry. 

We have studied the commensurability in I c (H a ) fur- 
ther by measuring a series of channel samples with dif- 
ferent values of the diamond cell length p and channel 
spacing s. Figure [3fa) shows I c (H a ) at T = 2.78 K 
for p — 0.5, 1, 2 /im, where all three sets of channels 
had s = 20 /im. For smaller p, the dominant peaks 
in I c shift to larger H a , although more complex oscil- 
lation patterns develop as well. Nonetheless, the Fourier 
transforms of the I c (H a ) data indicate that the lowest 
characteristic frequency in the spectrum for each p, ui° H , 
varies linearly with p [Fig. (3Jb)]. This provides evidence 
that the I c (H a ) peaks are indeed related to a match- 
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FIG. 4: Magnetic hysteresis in I c (H a ) for larger field sweeps, 
with H m ax as indicated for s — 20 /mi, p — 2 /on, T = 2.78 K. 
Curve for Hmax ~ 6.2 Oe is the same as in Fig. [2] with no 
hysteresis. 



ing of the vortex distribution to the constriction lat- 
tice. Because each vortex corresponds to one $o of flux 
($ = hc/2e rs 20.7 x I0~ 8 G - cm 2 ), the change in flux 
density in the channel AB c h that is required to add one 
vortex to each diamond cell will be determined by the 
area occupied by this flux. For widely separated channels 
(s > A^ bN ), the flux will extend ~ A^ bN into the banks 
on either side of the channel, while along the channel, 
the relevant length for the flux is p. Thus, one arrives 
at a rough estimate, AB c h ~ $o/2pA^ bN . However, if 
s is not large compared to A^ bN , the resulting overlap 
between vortices in adjacent channels will lead to an un- 
derestimate of AB^. In Fig- OH ) we see that a line 
through the uijjip) values has a slope that is approxi- 
mately 0.6 x 2A^ bN /<I>o, thus a somewhat larger AH a 
is required to achieve a particular AB c h- This is likely 
due in part to neglecting the overlap between vortices 
(s = 20 /im in this case), but is also related to the edge 
barrier mechanism. For a superconducting strip geome- 
try in a perpendicular field, B will be somewhat smaller 
than H a due to screening effects until H a 3> H s , where 
H s is the surface entry field [211 ]. 

In the opposite limit, s <C A r J bN , vortices in adjacent 
channels will be highly overlapping and the flux den- 
sity required for a one-vortex change becomes AB c h ~ 
&o/sp. We have varied the channel spacing s and ob- 
served the influence on I c (H a ), using s — 5, 20, 50 /im 
with p = 2 /im and T — 2.78 K for all three sets 
[Fig. OUc)]. The peak structure shifts to larger H a for 
smaller s, and the plot of oj° h vs. s in Fig. O^d) follows 
the trends described above, indicated by the dashed and 
dotted lines included in the plot. The s = 50 /im data ap- 
proaches the expected uj° h for widely separated channels, 
while the s — 5 /im data is close to the limit of highly 
overlapping vortices. In both cases, one expects a reduc- 
tion in (Jjj somewhat below AB~^ because of the edge 
barrier. A detailed calculation of the flux distribution in 
the channels, accounting for the channel structure, the 
two different superconductors, and the strip geometry, is 
beyond the scope of this paper. 
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FIG. 5: (Color online) (a) Example of dynamical hystere- 
sis in IV curve in the vicinity of one of the matching peaks 
(red/blue) compared to a reversible IV curve (green) [s = 
20 jLtm, p = 1 /im, T = 2.70 K]. (b) Density plot of the differ- 
ence of the flux-flow voltage between the outgoing and return 
current sweeps [AV = V ut(I) — V re turn(I)] as the color scale 
for different H a ; corresponding I c (H a ) superimposed (yellow). 



At the edge of a superconducting strip, vortices will 
enter when H a reaches H s , corresponding to the current 
density at the edge reaching a critical level, typically 
of the order of the Ginzburg-Landau depairing current 
density. Applying the standard edge barrier expression 
for H s [19j with our estimated film parameters leads to 
^NbGe _ 2.6 - 0.7 Oe, although the entry field into the 
ends of the NbGe channels is likely somewhat smaller 
than the Hf hGc estimate when one accounts for cur- 
rent distortions at the channel ends. Indeed, we typi- 
cally observe the first entry of vortices into the channels 
followed by oscillations in I c (H a ) for H a ~ 1 Oe. Per- 
forming a similar estimate for vortex entry into the NbN 
banks yields Hf hN ~ 8 Oe. We can probe the possi- 
bility of vortex entry into the NbN by increasing H a to 
progressively larger values H max before reducing it and 
checking the reversibility of I c (H a ), as vortices trapped in 
the strong-pinning NbN will exhibit an irreversible mag- 
netic response and will offset the net magnetic field ex- 
perienced by the vortices confined to the channels. For 
small H max , I c (H a ) retraces completely [Fig. Ufa)], cor- 
responding to the entry of vortices only into the NbGe 
channels. However, for H max > 8 Oe, I c (H a ) becomes 
hysteretic, with the opening of the hysteresis loop grow- 
ing with H max (Fig. 0]). Also, the matching peak struc- 
ture on the return branches of I c (H a ) becomes washed 
out for larger H max , as the disordered distribution of vor- 
tices that occurs in the strong-pinning NbN when H a is 
reduced randomizes the potential for the vortices moving 
in the channels. 

For H a below the threshold to introduce vortices into 
the NbN banks, in the vicinity of the I c {H a ) matching 
peaks, we often observe a completely different type of 
dynamical irreversibility consisting of hysteresis in the 
IV curves. Figure |5J a) shows two example IVCs for the 
outgoing and return current sweeps, one between match- 



ing peaks with no hysteresis, the other near the second 
matching peak with clear hysteresis. Thus, in this second 
example, there is a clear irreversibility in the flux-flow 
voltage, but not in the critical current. 

We can combine all of the measured IVCs for a par- 
ticular channel configuration and T by making a density 
plot, where the color scale is the difference between the 
voltage on the outgoing and return current sweeps. We 
superimpose the corresponding I c {H a ) curve for reference 
[Figure [5[b)] . This particular example, with s = 20 /im, 
p = 1 /im, T = 2.70 K, shows regions of dynamical hys- 
teresis near the first three I c (H a ) peaks. Over the range 
of drive frequency that we have studied, 20 — 400 Hz, we 
observe no change in this response. The upper limit of 
the data on the current axis is set by the point where 
the flux-flow voltage approaches the Larkin-Ovchinnikov 
instability point [22j , where the channels switch abruptly 
to the normal state. 

This hysteresis in the IVCs may correspond to a dis- 
tortion of the vortex distribution as the driving current is 
reduced that allows the vortices to keep flowing at higher 
velocities than when the current was initially increased. 
We note that not all of the diamond channel measure- 
ments displayed this dynamical hysteresis. In particular, 
none of the channel configurations that we have studied 
exhibited this type of irreversibility when measured at 
the highest temperature of our experiments, T = 2.90 
K. Figure [B] contains example IVCs and a similar density 
plot with superimposed I c (H a ) curve to that in Fig. Efb) 
for a set of channels with s — 50 /im, p = 2 /im, measured 
at T — 2.90 K where there is no evidence of dynamical 
hysteresis. This may be due to the change in the inter- 
vortex interaction strength as T approaches T^ hGc . We 
are currently investigating this dynamical hysteresis in 
our channels further. We note that hysteretic dynamics 
for vortices in periodic arrays of antidots were recently 
reported [23[. These were connected to previous theo- 
retical work involving the transition to turbulent flow 
related to the interplay between interstitial vortices and 
those pinned in the antidots [HI, H3] ■ The origin of the 
hysteresis in our system is likely somewhat different, as 
all of the vortices are confined to the weak-pinning chan- 
nels. 

In summary, we have measured vortex dynamics in 
weak-pinning channels containing periodic constrictions 
that are small compared to the vortex size. Over much 
of the magnetic field range that we have studied, all of 
the vortices are confined to the channels and the chan- 
nel structure results in strong matching effects between 
the vortex distribution and the constriction lattice. In 
the vicinity of the matching peaks, we often observe a 
dynamical hysteresis in the vortex response that may be 
related to a distortion of the vortex distribution. 
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by the National Science Foundation under Grant DMR- 



5 




0.00 0.10 0.20 0.30 2 4 6 8 10 12 
I(mA) H a (Oe) 

FIG. 6: (Color online) Example of measurement at higher 
temperature with no dynamical hysteresis present [s = 50 fim, 
p = 2 fim, T — 2.90 K]: (a) sample IV curves for outgo- 
ing and return sweeps together, one in vicinity of one of the 
matching peaks (0.68 Oe) and one in between matching peaks 
(1.6 Oe). (b) Density plot for this channel configuration and 
temperature calculated with same technique as in Fig. (5jb); 
corresponding I c (H a ) superimposed (yellow). 
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